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Abstract 

This paper proposes a one-dimensional lattice model with long-range interac- 
tions which, in the continuum, keeps its nonlocal behaviour. In fact, the long- 
time evolution of the localized waves is governed by an asymptotic equation of the 
Benjamin-Ono type and allows the explicit construction of moving kinks on the lat- 
tice. The long-range particle interaction coefficients on the lattice are determined by 
the Benjamin-Ono equation. 

1 Introduction 

This paper is concerned with the effects of long range interactions in the behaviour of 
solitons (kinks) in a lattice system. In particular, we show that the introduced lattice 
model asymptotically leads to a nonlocal continuum equation. By nonlocal we mean the 
model where the localized solution is of order of the scale parameter [it parameter which 
has the dimension of length), and thus it is physically acceptable to consider wavelengths 
comparable with the scale parameter. Its equation of motion contains integral, integro- 
differential or finite difference operators in the spatial variables, while the wave propagation 
velocity depends on the wavelength. In the present model there are two characteristic 
lengths: (a) the lattice spacing and (b) the radius of particle interaction range. 

The nonlocal models can be divided into two classes: discrete and continuous. For 
specific cases, it is simpler and easier to use the quasicontinuum description of the discrete 
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medium. Its essence is an interpolation of functions of discrete argument by a special class 
of analytic functions in such a way that a correspondence condition between the quasicon- 
tinuum and the discrete medium is fulfilled. The advantages of such an approach is the 
description of discrete and continuum media within a unified formalism and their correct 
generalization. In that case the quasicontinuum model is also applicable to macrosys- 
tems. Analytic solutions valid for wide nonlinear excitations have been obtained using 
the continuum approximation in [[[], |J and the quasicontinuum approximation in || f|, f|. 
In particular in [|, it has been shown that very narrow soliton-like excitations can 
propagate without any energy loss due to discreteness effects. 

The lattice model we study, possesses an inherent nonlocal behaviour, ie the motion of a 
particle at node n depends on the motion of its neighbours. Such models for different kinds 
of particle interactions have been studied in a series of papers J7], So far the classical 

procedure for studying them, at the continuum limit, is to assume slowly varying fields 
and expand them using Taylor's series. Then, approximate partial differential equations 
for the continuous fields were obtained. 

In this paper, we show how to keep the nonlocal behaviour of the lattice model at the 
continuum level, or in other words, how to describe the nonlocality at the continuum level. 
This transition from the continuum to the discrete system with long-range interactions 
and vice-versa, is straightforward due to the quasicontinuum approximation. The idea of 
preserving the nonlocality at the continuum limit occurs in condensed-matter physics to 
describe phase transition in crystals, nonlinear waves in crystals and biological molecules. 
For example, modified versions of the discrete nonlinear Schrodinger equation which keep 
the nonlocality at the continuum and describe problems in biomolecules have been studied 
in the past Q 0. 

It is worthwhile mentioning some interesting studies of anharmonic atomic chain includ- 
ing long-range interactions due to large Coulomb coupling between particles ]12|, [13], [14], [L5 



The long-range potential is of the inverse power-type or Kac-Baker form. In the case of 
inverse power type of the fourth-order, the continuum limit yields an integro-differential 



equation involving a Hilbert transform [y|. Using a perturbative technique the authors 



derive a mixed modified Korteweg-de-Vries with Benjamin-Ono equations ||15|| , which is 
close to the result presented in the following. 

The paper is organised as follows: In Section 2 we introduce the lattice model; while 
in Section 3, we deal with its continuum approximation and show that it leads to the 
Benjamin-Ono equation for long-time evolution. In Section 4, we derive the type of long- 
range particle interactions on the lattice using the fact that the lattice model asymptotically 
leads to the Benjamin-Ono equation. Finally, numerical simulations of the discrete model 
based on the Benjamin-Ono soliton are given and discussed in Section 5. 
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Figure 1: One-dimensional atomic chain of masses and springs with first and second neigh- 
bour interactions. 

2 The Lattice Model 



We consider a one-dimensional lattice with mass M, lattice spacing b and long-range in- 
teractions. The Lagrangian of the system is 

C = £ \mu\ - £ h{n - n') (u n - u n ,f - £ \k 3 (u n+1 - u n f . (2.1) 

n 1 nn' 4 " d 

Figure 1 presents a one-dimensional chain of masses and springs with second neighbour 
interactions. 

The equations of motion obtained from the variation of ( ]2.1|) are 



M 



MlL n = Y^i 771 ) ( U n+m + Un-m ~ 2m„) + k 3 («n+l ~ U n ) 2 - (u n ~ M n _i 



(2.2) 



m=l 



where k 3 is a positive constant (k 3 < 1). ^(n) are the interaction coefficients of the parti- 
cles, related to the force constants $(n) (see equation( [2.4|) ) which determine the effective 
characteristics of the elastic bonds between the particles. The lattice parameters \l/(n) can 
be obtained from the interaction potential between particles or atomic bonds in crystals 
(same as for $(n)); e.g. the Born-Meyer potential, the Morse potential, the Lennard- Jones 
potential, the covalent potential, etc. 

Here, we consider the chain to be homogeneous and so, the unique definition of $(n) 
leads to the identity $(n) = $(— n). If $(n) differs from zero for \n\ < A4, then every 
particle interacts with M. neighbours to the right and to the left. The case M. = 1 
corresponds to the simplest model of interaction between first nearest neighbours; while in 
real mechanical systems the action at a distance is always restricted, ie M. is finite. 

The distances between the particles do not change and therefore the forces acting on 



them are equal to zero [0. This leads to the following condition for the force constants 



$ (o) = -E> 



n 



(2.3) 
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where J2' n m eans summation over all n / 0. This equation is deduced from the lattice 
energy invariance; ie the invariance by translation of the chain as a whole. 

Assume now, that the particles interact nonlocally with measure ^(n). Then it can be 



shown that 16 



$(n)=-tf(7i), n^O; $(0)=^*(n). (2.4) 

Note that for real mechanical systems all \l/(n) > and hence $(0) > 0, $(n) <0 (n^ 0). 
However, a system may be constructed, for which some particles would have negative 
interaction coefficients; which is the case of our model described below. 



3 The Continuum Approximation 

The set of nonlinear differential equations for a discrete system is generally complex and, 
quite often, its continuum approximation is considered (which is the case here). By as- 
suming that the displacements are slowly varying over the lattice spacing, the continuum 
approximation may be derived as follows: 



3.1 The Nonlocal Model 



The Lagrangian ( 2T[) in the form which is invariant with respect to n and x representations, 
can be written as 

C = - < pu,ii> +]- < u, $m > ~ < Au, k 3 Au 2 >, (3.5) 

Zi Zi -_> 

where p is the mass density, while u = u n (t), Au = b~ l {u n+ i — u n ) and u = u(x, t), Au = u x 
in the discrete and continuum space. 

The various terms in ( p.5| ) for both the discrete and the continuum representations are 
given as follows 

nn' 

— — ( $(x — x') u(x) u(x') dxdx', (3.6) 

- < Au, k 3 Au 2 > = i h{u n +i ~ u n ) 3 
6 6 n 

= i / huldx. (3.7) 
3 Jr 

The above forms contain elastic (long-range) interactions for both the discrete model and 
its quasicontinuum counterparts. 
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The equation of motion of the medium in the (x, t) representation is 



pii(x,t) = — <£>(x — x')u(x', t) dx' + k^(u x ) x . (3 



The corresponding condition (|2.3| ) in the x representation becomes J R $(a; — x') dx' = 
which implies that $(x — x') can be represented as 

$(z - x') = tp(x) 5{x - x) - V(x - x% (3.9) 

with ty(x — x 1 ) = ty(x' — x) and ip(x) = J ty(x — x') dx' . 
Its quasicontinuum counterparts are of the form 

$( n - n ') = ip(n) 5(n - n) - *(n - n'), (3.10) 

with $?(n — n') = ty(n' — n) and ip(n) = J2 n ' ~ n ')- Thus ^>(n — n') and ^>(x — x') are 
the interaction coefficients and the interaction coefficient density of the model connecting 
the points (n, n') and (x, x'), respectively. 



In the x representation the force constant can be expressed as JTG 



^x-x') = ^c(x-x'), (3.11) 

where c(x — x') is the so-called kernel operator of elastic moduli. 
Then, equation ( |3.8p after integration by parts becomes 

pu(x,t) = ^ J c(x - x')^^- dx' + k 3 (u 2 x ) x , (3.12) 

which is quite similar to the equation obtained from the nonlocal elasticity of the one- 
dimensional linear case (ie with the k$ term vanishing) [T7], TB 



In this case, by using the kinematical definition for the strain e(x',t) = du(x',t)/dx' 
and the nonlocal constitutive equation for the stress a(x, t) = f_ t c(x — x') e(x', t) dx' where 
I is the interaction range length, equation fl3.12| ) can be recast into 

pu(x,t) = a x , (3.13) 

which is the equation of motion of a one-dimensional elastic medium. [Note that, the 
integral over a finite radius of interaction defines the finiteness of action at a distance 
determined by this radius (that is, I = Mb). Specifically, $(y) is not zero for y G [—1,1] 
and it equals to zero otherwise; while, if the chain is bounded within an interval 2L then 
/ « L] 

In this connection, equation ( |3.13| ) is the starting point in continuum mechanics and 
expresses the balance of linear momentum. It is endowed with appropriate constitutive 
equations for the stress leading to partial differential equations for the determination of 
the displacement field. For example, a simple linear elastic model of the form a = Ee with 
E denoting the Young modulus (identified as Cq below) leads to the classical wave equation 
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u tt = c 2 u xx (c 2 = E/p). If gradient-dependent nonlinear elasticity model is adopted |H5 
of the form a = /(e) — ce xx the appropriate equation of motion reads f{u x )u xx — cu xxxx = 
pu tt which contain the k^{u 2 x ) x term listed above, even though a more interesting case 
results when f'(u x ) is non-monotonic. Other types of such constitutive models and the 
corresponding equations of motion are discussed in ||20|| . A case of particular interest in 



terms of its discrete model counterpart arises when viscosity is considered. 

For a simple Kelvin- Voigt model of the form a = Ee + pe where p is the viscosity and 
e the strain rate, the resultant equation of motion with A = p/p is u u = c 2 u xx + Xu xx , ie 
the damped wave equation. It is an open question to examine a Lagrangian-like lattice 
version of this model. 

On the other hand, by letting 

c = c 2 5(x-x') + ci(x-x'), (3.14) 
where c corresponds to the local response with Cq denoting the local elastic modulus and 



Ci describing the nonlocal behaviour of the model, equation (|3.12|) simplifies to 



pu(x,t) = c 2 u xx + ^- J r Ci(x - x') du ^ f - dx' + k 3 (u 2 x ) XJ (3.15) 



which is further examined below. 



3.2 The Long-Time Evolution of the Localized Wave 

In this section we investigate the asymptotic behaviour of equation ( |3.15| ) for large times. 
By assuming that the contribution of the nonlinear term is significant throughout we rescale 
the nonlinear and the dispersive term by introducing a small parameter e < 1, as follows 

2£; 3 = 7£, c 1 = eg. (3.16) 

Here, 7 is the rescaled nonlinearity coefficient and g is the rescaled elastic kernel. Then, 
for large times, the asymptotic expansion for the displacement field can be expressed as 

u(x,t) = u (^t) +eu 1 (x,t) +0(e 2 ), (3.17) 

where £ = x — v t is a shifted coordinate and r = et is the slow time variable. 

Using equations ( |3.16| ) and ( |3.17| ) and keeping terms of order e only, equation fl3.15|) 
becomes 

r 

2pvu 0(T + -g^J R 9(£ - O u o(> d£' + 7 u 0li uo^ = -c 2 u lxx + p u ltt , (3.18) 
where Cq = pv 2 . 

Then, the secularity condition implies that both the left and right-hand sides of 
( [3. 18] ) are zero. Thus, the right-hand side leads to the standard linear wave equation for 
U\\ while the left-hand side gives the long-time behaviour of uq, ie 

d f 

2pv Uo^t + ^7 J R 9(£~ O M oe d£' + 7 m « u Q ^ = 0. (3. 19) 
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Finally, by letting 

iK* = .F, g(z) = ^G(z) } (3.20) 



where z = £ — £' and y = £ — z, equation (|3.19|) becomes 

2 P vF T + 1 FF l : + J R G(Z-y)^-F(y,T)dy = : (3.21) 



which transforms to the Benjamin-Ono equation for 



X -T, ^=il7, G(£-y) = ---!-. (3.22) 



2pv 7 ' 7r £ — y 

Note that, under the swap of the role played by x and t in (|3.21|) , we arrive at an 
integro-differential equation usually met in nonlinear wave propagation in active media 
with dissipative process p2| . The latter is due to viscoelasticity described by an integral 
law involving exponential kernel in the integral part of the equation of motion. This kernel 
is then a function of time introduced by a relaxation process. 

3.3 The Benjamin-Ono Equation 

The Benjamin-Ono evolution equation 

U T + UU i + aH(U £i ) = 0, (3.23) 

with a being a positive parameter, was originally derived in [^3|, ^] for interval-wave 
propagation in a two-layer system: one shallow and the other infinitely deep. In this 
context, {/(£, T) represents the amplitude of the interfacial wave produced by an initial 
disturbance, say {/(£, 0). The operator H designates the one-dimensional Hilbert trans- 
form, ie 

a r 1 



a H(Ux) = - / Uyy dy, (3.24) 

7r Jr y — 4 

which implies that ( |3.23|) is an integro-differential weakly nonlinear evolution equation. 

Equation ( |3.23| ) has a simple solitary wave solution in the form of a Lorentzian (alge- 
braic) shape pi], p6| , ie 

where Av is the amplitude of the wave and ^ measures the wavelength. The velocity of the 
soliton is amplitude-dependent and equals to v. 

Remark: Regarding the long-wavelength limit of the lattice model, the classical con- 
tinuum approximation (ie Taylor's expansion for the fields) of (|3.15 ) leads to the local 
nonlinear partial derivative equation of the Boussinesq type |2j 

lift ~~ (? u xx + dUxUxx + bu xxxx = 0, (3.26) 
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where the parameters c, a and b depend on the atomic interactions. Alternatively, equation 
( p.26| ) may be directly deduced as a special case from the continuum approximation of (|2.2|) . 
The long-time evolution of ( 3.26|) is given by the Korteweg-de-Vries equation, 



V t + VV X + V xxx = 0, (3.27) 

when multi-scale techniques (or perturbation methods) as described in Section 3.2 are used. 
In addition, the Benjamin-Ono equation ( |3.23| ) can be transformed into the Korteweg-de- 
Vries equation by replacing aH{U^) by however, this transition is not an obvious 
one. 



4 Evaluation of the Interaction Coefficients 

In this section we derive the lattice interaction coefficients \&(n) of the model ( [2.1|) using the 
fact that its asymptotic expansion, for large times, leads to the Benjamin-Ono equation; 
and therefore, the kernel of the operator of elastic moduli c(x) is known. 

Since there is a correspondence between functions of discrete argument and a certain 
class of analytic functions, the interaction coefficients can be evaluated either in the space 
of functions of discrete argument or in the space of analytic functions. However, since all 
these spaces are isomorphic to each other the interaction coefficients are obtained using 
the x representation (due to simplicity). 



Equation ( |3.11| ) by setting x — x' —> x becomes 



d 2 

0; 

■2 



= -cl5"(x)-eg"(x), (4.28) 

where ( |3.14| ) and ( |3.16| ) were used. But from the long-time evolution analysis of the 
localized wave the explicit form of the function g is known. Thus, using ( |3.20|) and ( |3.22|) 

we get 



$( x ) = - c l5"{x)-eG"'{x) 

6ea 1 



Due to Fourier transform, the force constants in the space of discrete argument are 



related with the ones in the space of analytic functions as |T6f: 



$(n) = $(n&) = / $(x) $ b (x - nb) dx 

JR 



R 

cll 1 (n)-—l 2 (n), (4.30) 

7T 
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where 8b(x) = (2tt) 1 J B e tkx dx = sin(7rx)(7rx) 1 and 

TAn) = 



sin[f(x-7i6)] 

o (x) 7 7- — ax, 

/>• 7r(x — no) 



T 2 {n) 



1 sin[|(x — nb)] 



dx. 



Irx a ir(x — nb) 

Thus, it is a matter of algebra to evaluate the two integrals and obtain 



(4.31) 



1,(0) 

Mo) 



Ti- 



ll 



3 6 3 ' 

4 



24 6 4 ' 



(4.32) 



while for n ^ 



Zi(n) 



2 (-l) n 
6 3 n 2 

7T 2 (-1)" " 

¥ n 2 



+ 



•l) n 



n 2 7T 2 



(4.33) 



Therefore, the lattice force constants are given by 



2^2 



$(0) 
$(n) 



36 3 
2cg( 



-l) n 6aen (-1)' 



+ 



7T, 2 7T 2 



n ^ 0. 



(4.34) 



b 3 n 2 6 4 n 2 

Finally, the lattice interaction coefficients \P(n) in (|2.2|) (due to (|3.10|) ) are given by 

*(n) = $(0)<J(n) - $(n). (4.35) 

As one can observe from Figure 2, by keeping the nonlocality of the discrete model to 
the continuum one, and specifically, the choice of an integrable nonlocal equation as its 
asymptotic limit for large times, leads to the construction of a lattice model with negative 
interaction coefficients between its particles. In fact, Figure 2 illustrates the change of sign 
of ty(n) and its variation with n for an atomic chain where each particle interacts with its 
ten neighbour particles. [The values of the other parameters, ie (a, b, e, Cq) are listed in the 
figure caption.] 



5 Numerical Simulations 

The behaviour of genuinely discrete systems is quite different than the one of their contin- 
uum counterparts and these differences have been studied in numerous papers both from 
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Figure 2: The function *(n) given by ( ggg )-( gg5p for a = 1, v = 1/4, e = 0.1, n = 125 



and phase velocity (a) c = 0.6 solid line (b) c = 0.65 dashed line (c) c = 0.7 dotted line. 

the mathematical and the physical point of view. Many of the methods used and their 
physical implementation in models with nearest neighbour interactions can be found in 
p7| , |28| , |29j. In addition, in |30], |31| the effects of discreteness in solitary waves in nearest 
neighbours lattices have been mathematically investigated. 

In our case, the lattice model keeps its nonlocal behaviour at the continuum when its 
asymptotical limit for large times is the Benjamin-Ono equation. However, it is not obvious 
if the solution of the Benjamin-Ono equation will be the solution of the lattice model, as 
well. Thus, in this section we investigate the dynamical behaviour of the lattice model 
( [2. 2D numerically, using as an initial condition the soliton solution given in Section 3.3. 

We propose a numerical scheme by directly considering the lattice equations ( |2.2| ), 
employing a Runge-Kutta method of fourth order and imposing pseudo-periodic boundary 
conditions, ie 

u(N + i) = u{i) + M , ie[-M,M], (5.36) 

where N is the number of particles in the lattice, M. is the range of particle interactions 
considered and uq is the amplitude of the initial kink. 

The initial conditions for the displacements and the velocities of each lattice particle 
are given by the analytic expression 



lea 

u[x,t) = — — arctan 




(5.37) 



and its time derivative, respectively. Here xq denotes the soliton position and we have 
set the mass density equal to one (p = 1). This solution represents a kink with thickness 
A = anr/u, ie the number of particles participating in the kink solution is proportional to 
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(b) 



Figure 3: (a) Initial soliton solution (|3.25|) and (b) Initial kink solution ( |5.37| ) for a 
v — 1/4 and n = 125. 



the kernel factor a and decreases as the wave parameter v increases. [The solution (|5.37 ) 
has been obtained from the solitary wave solution of the Benjamin-Ono equation fl3.25|) 
after using the change of variables and fields given in Section 3.2]. Figure 3 illustrates (a) 
the soliton solution U(Q with ( = £ — vT given by ( |3.25| ) and (b) the kink solution u[x) 
given by (|5.37| ) for the choice 6 = 1. 

In all our simulations we have chosen the following fixed values for: k 3 = 0.1, e = 0.1, 
6=l,a = l,M = l. However, different values for the above parameters give qualitatively 
the same results. The initial position of the kink is chosen equal to n = 125, while the 
value of the kink velocity cq is chosen so that the first interaction coefficient is positive (ie 
\l/(l) > 0) and the total number of particles in the lattice is equal to 250. 

Figure 4 summarises the main results of our numerical simulations when models with 
different number of non-nearest neighbour interactions are considered and different values 
for the velocity and wavelength of the initial kink have been chosen. In particular, Figures 
4(a), 4(e) and 4(f) present a solution corresponding to a lattice model with first, tenth 
and twentieth neighbour particle interactions when the initial kink wavelength and phase 
velocity are: v = 1/8, c = 0.7; v = 1/4, c = 0.6; and v = 1/8, c = 0.6, respectively. 
In each model, the kink contains from twelve up to twenty four particles, which means 
that the continuum approximation is valid. [Recall that the kink thickness depends on 
the parameter z/, due to fl5.37p ]. On the other hand, Figures 4(b), 4(c) and 4(d) present 
the kink dynamics for lattice models with second, fourth and seventh neighbour particle 
interactions. However, we choose the kink width in all three cases to be small, ie the kink 
contains six particles with v = 1/2, cq = 0.54 in the first case, and three particles with 
v = 1; Co = 0.6, Co = 0.7 in the last two. 

From Figure 4, one can observe that the initial soliton solution of the Benjamin-Ono 
equation is also a solution of the lattice model with long-range interactions given by (|2.2| ). 
In fact, the initial lattice kink relaxes and propagates in the lattice space with small 
oscillations. It can be observed from Figures 4(a), 4(e) and 4(f) that as the kink thickness 
increases (ie v decreases) the lattice radiation is eliminated. Numerical simulations of the 
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discrete nonlocal model for kink widths equal to three or six lattice spacings have been 
presented in Figures 4(b)-4(d): at this level the kink changes slightly and loses energy 
by radiation of small amplitudes waves. In all cases where the kink contains at least 20 
particles, it is remarkable stable and emits (almost) no radiation. 

Let us conclude by saying that, in the long- wavelength approximation (ie, when the 
average kink wavelength is larger compared to the lattice spacing) the initial (continuum) 
kink, which has been derived using the quasi- continuum approximation of the discrete 
model, relaxes to a stable discrete kink without any oscillations and our method is accurate. 

6 Concluding Remarks 

The present work is concerned with the dynamics of lattice models with long-range and 
nonlinear interactions. The main objective is the description of the dynamics of the lattice 
model, at the continuum level, by keeping the nonlocal nature of the discrete system. Our 
model is the simplest one consisting of a one-dimensional chain of point masses connected 
by springs. We have developed a theory for such a discrete system and introduced the 
notion of quasicontinuum which allows us to treat the discrete and continuum models in 
the framework of the same formalism. Attention has been focused to the analysis of an 
approximate model and its transition to nonlocal elasticity in the limit of long waves. In 
this formalism, we have set up a one-to-one correspondence between functions of discrete 
arguments and a class of analytic functions. The procedure is also applicable for integral 
operations on these functions. The method brings out a resemblance to the problem 
of interpolation of a discrete function by a smooth function satisfying some smoothness 
conditions. The technique allows us (Section 3) to use the same representation of the 
Lagrangian for all cases: discrete and analytic. 

The equation of motion of the model in the continuum representation is similar to that 
of a nonlocal model of elasticity including a local nonlinear term. Nevertheless, such an 
equation of motion is not easy to deal with and a perturbation method has been employed. 
More precisely, we have considered the long-time evolution of nonlinear waves traveling in 
the lattice. The asymptotic method has been used to transform the equation of motion 
into an equation of the Benjamin-Ono type which governs the long-time evolution of the 
initial signal. For a particular choice of the nonlocal elastic kernel or interaction force 
function, the dispersive operator of the equation is a Hilbert transform corresponding 
to the Benjamin-Ono equation. In this particular situation, a localized wave solution of 
soliton type has been obtained. In order to test the analytical conjectures thus developed, 
some numerical simulations have been performed directly on the lattice model for different 
ranges of nonlocal action. It has been observed from numerics that the initial soliton 
(kink) solution propagates almost without perturbation and it is very stable. Moreover, 
wide kinks propagate in the lattice without any loss due to discreteness radiations. 

Such a study can be extended to other kinds of nonlinear terms (double-well or sine- 
potential, higher-order interactions) or to arbitrary particle interactions. Extension to two- 
dimensional lattices can also be envisaged. The problems of kink collision or the influence 
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(a) 



(b) 




Figure 4: Kink plane for (a) one-range interaction, v = 1/8, c = 0.7 (b) two-range 
interaction, v — 1/2, Co = 0.54 (c) four-range interaction, v — 1, cq — 0.6 (d) seven-range 
interaction, v = 1, c = 0.7 (e) ten-range interaction, v — 1/4, c = 0.6 (f) twenty-range 
interaction, z/ = 1/8, c = 0.6. 
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of perturbations (forces, damping, lattice defects) on the kink motion is of interest. Some 
of these problems will be addressed in the future. 
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